UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER LABORATORY 

NEW COMPUTER LIBRARY ROUTINE B3-1XP1-25 



TETLE: 

TYPE: 

LENGTH: 



exponential 

closed; relocatable T mnemonic 

30 words 



TEMPORARY STORAGE: ten words at fixed memory locations 0-9 and the last two 

words of the program « 

approximately 800 microseconds (October I96?) 



D10RATION: 

FAST REGISTERS 
CHANGED: 

PARAMETERS: 

ACCURACY: 



BISCRIFTION: 



COMPARE: 



none 

link in Ml? 



^-kO 



relative error < 2 for -ink < % < mk 

This subroutine replaces the number z in the accumulator by 

if z < -6k &nk < -87 
e 2 if -6k Ank < z < 6k itrk 

and sets OV if z > 6k £nk > 87 

For < x < H>vk the formula e x = 1 + x|l-x/2+F) is used 

where F is the continued fraction F = r^ + fi + Fi "*" ""9 

a «* x f(k{kn - l)), which is computed to eight terms . 

For a > Ankg e » k e , where h = integer part {zf&nk};, 
and < x < take For negative values of z, e a is computed 
as l/e" , 

No Macomj, "On the computation of exponential and hyperbolic 
functions using continued fractions „" JACMj, Oct '55 



BATE: November lk f 1962 
PROGRAMMED BY: J. Nievergelt 



.1- 



EXP 



I 



SFR7*0 

CAM13*0 

CMM15«S,1 



4 



5 



-<JOC2,0,2R 



(una,* 

»STHlg,0 



1 



>S?R1S,X 
DlV8 t 3 f m 

Sy89»3 g 64 



-JDC3,1,5R 






t 



DtV15,3 
-JDCO, 2,8511 



AD09,3 9 64 

SOAlfoO 



€ 



10 



12 



AN08,3,26R 
HPY8,3,I7R 



STRXS,! 
MPYX0,3,IQ48 

STRIS^O 



STR8,3^18R 
DIV9»3,3 



STR8 e 3,f9R 
CAD9 a 3,4~ 



SYR15,1 
CAD9*3 fl 3^ 



STRlSjl 
CA09,3,1> 

STRisq 



STRS^S^ 

CAD§S«3> 

STRaS,l\ 



« address of Z 

jump If Z > 

if Z negative, (MIS) « 1 
and replace Z by j ZJ 

j Z | <o memory 1 
Z/fn4 



jump If Z/ln4 < 64 

A<seu « e* » 

jump If Z was negative 

set OV If Z was positive 



reset Z/ln4 

h » Integer part (Z/in4) « Mil 

fractional part f In eccu 

x a f Sn 4 
x -«> memory S 

x/2 <© memory 3 

x^/4 «* memory 2BR 



a * (x /4|/3 -0 memory 29ft 
initial value ^""s 4 * memory 4 
Initial value d e 3 o memory 5 

initial value 

A-l * Bq • ^ 

Initial value 
Aq - B -3L - 



©- 



35 



16 



n 



18 



19 



80 



21 



m 



is 



24 



ffl-^r 



86 



587 



511115,0 
CSMl4,f,?/ 



-£~ CAM25,a,6 



CAD8,3,29K 



MPVlS.l 
"Xt)Dl5,6 
XCH15,1 
STR8,5,6 



CA08,3,29R 

MPY15,1 
ADDlSjJ) 



XCH15,0 
STR8,3,8 

CA09,3,8 



CAM25,f,4 

ADP15.0 

STR15,1 
ADDIS ,0 
STRl5 f l 

VliD8,5,S8ft 



STR8,3,29R 
CJUl4,S,i3R 



CA08,3,T 
DI1V8,3,9 



ADP9,3,1 
SUB8,3,3 



VI 08,5,3 

ADD9,3,; 



ADE12,0 

JZM13,2,g5R 



V8D9,3,I 



yiFR7,0 
JLH15.0 



OCT© ,17777 ,Vmi ,17777 



loop counter 

(MIS) - address of A R ^g 

a 



n 



B A n=t + A n~i * A 
old A n „i - A n ^ 

a 



a *iHe + B n~l - A n 

old B n ~i - B n -3 



a" * 8 

(MIS) » address of a' 

A 1 ' <-» A* + 8 

d <=— d ■$• a' 

(x2/4)/d - a 



(M15] 



F ■ A/B 

F + 1 



(x/t) 



+ x/(F + 1 - x/S) 
jump Iff Z was positive 



e* s 4,n e 



0Ct02S13, 110t? , X7372 ,03401 



!ra 4 



A n«t < 



